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Abstract 

To leading order approximation, the physical content of generalized parton distributions (GPDs) 
that is accessible in deep virtual electroproduction of photons or mesons is contained in their value 
on the cross-over trajectory. This trajectory separates the t-channel and s-channel dominated GPD 
regions. The underlying Lorentz covariance implies a correspondence between these two regions 
through their relation to GPDs on the cross-over trajectory. This point of view leads to a family 
of GPD sum rules which are a quark analogue of finite energy sum rules and it guides us to a new 
phenomenological GPD concept. As an example, we discuss the constraints from the JLab/Hall 
A data on the dominant w-quark GPD H. The question arises whether GPDs are governed by 
some kind of holographic principle. 
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1 Introduction 



Generalized parton distributions (GPDs) [H [21 [3] are partonic amplitudes which contain the non- 
perturbative physics of perturbatively factorized exclusive leptoproduction amplitudes [HE]. De- 
pending on the parton momentum fractions, a GPD is interpreted either as the probability ampli- 
tude of meson-like t-channel exchange or as s-channel exchange of a parton. Compared to parton 
distribution functions (PDFs) and distribution amplitudes, they reveal the partonic content of 
hadrons, protons in particular, from a complementary perspective; see for reviews [SI EJ- For 
instance, the partonic angular momentum, naturally defined by the decomposition of the angular 
momentum tensor in terms of quark and gluonic degrees of freedom, is given by the first Mellin 
moment of certain GPDs [8]. Quantification of such partonic decomposition of the proton spin is 
one of the main goals for measurements of GPD related processes. 

During the previous decade, experimental effort has been devoted to the measurement of deep 
virtual Compton scattering (DVCS) |3 dm dH [T21 H31 H31 ESI HSl IIZl CHI HSl I2D] axid production 
of mesons, see e.g. Refs. [2~Tj 1221 [231 I2H [25], yielding an increasing amount and precision of 
experimental data. On the other hand, GPDs are intricate functions, depending on the parton's 
longitudinal momentum fraction (x), the skewness (rj), the momentum transfer squared (t), and 
the factorization scale (/x 2 — usually set equal to the virtuality Q 2 of the probe). This makes 
their extraction from experimental measurements, where they enter as a convolution with a hard- 
scattering amplitude, a rather awkward task. Since the partonic momentum fraction dependence 
of GPDs is integrated out in the amplitude, the degrees of freedom that can be constrained by 
experiment do not provide the momentum fraction shape of GPDs for fixed 77. Consequently, GPD 
moments, including the angular momentum of quarks, cannot be directly revealed. 

The analysis of deep virtual exclusive processes is based on GPD ansatze and models which take 
into account that GPDs reduce by sum rules to form factors and in the forward limit to PDFs. As 
the basic GPD ansatz, written in essence as a product of form factor and PDF, the popular VGG 
code [26\ successfully describes the first DVCS beam spin asymmetry measurements at HERMES 
[TO] and CLAS [XT] , but it fails in the case of more recent data on the unpolarized DVCS cross 
section, measured by the Hall A collaboration [T3]. Other GPD models [27J [28j [29], [30] are similarly 
(un) successful in describing DVCS observables measured in fixed target experiments. In collider 
kinematics, no published GPD model, which is (a) mathematically consistent, (b) has a reasonable 
t-dependence, and (c) satisfies the evolution equation, is able to describe the DVCS data at leading 
order (LO) accuracy. Taking into account radiative corrections, the measured DVCS cross section 
can be described, as advocated in Refs. [3TJ [32], and fitted in Refs. [331 [31]. Also, Regge inspired 
modelling of off-shell amplitudes [35] can be tuned to be consistent with experimental data [3"6ll3~T] . 
If unitarity constraints are taken into account in the fixed target kinematics, the model [3S] is in 
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fair agreement with the unpolarized DVCS cross section measurement [38]. Thereby, it provides 
a physical interpretation in terms of hadronic degrees of freedom, where the desired view inside 
the nucleon cannot be delivered. 

Thus, although GPD modelling is guided by the relation of GPDs to PDFs and form factors, 
by Regge phenomenology, various model calculations, lattice simulations etc., the resulting models 
are only partially able to describe experimental data. This, in combination with the wide-spread 
implicit belief that the resulting GPD models might be realistic in a more general way, leads then 
to endless discussions, such as 'Which model is right, A or B?' (when both are maybe wrong) or 
'GPDs or not GPDs?', to claims that GPDs can be constrained from present lattice data, and to 'a 
conjectured proof of the breakdown of the GPD formalism (for some details see [3H]). Sometimes, 
criticism of GPD models is unjustly extended to GPD representations. In our opinion it is unclear 
whether sufficient understanding of GPDs is reached in the proposed models. Below we shall 
demonstrate that there should be no problem to describe presently measured DVCS observables 
within the GPD formalism. Hence, we should question the present phenomenological approach, 
which is based on GPD models that we consider to be ad hoc constructions. 

We would like to come back to two basic questions, which should be clearly answered. 

• What GPD information can be revealed from experimental data? 

• What do we learn from this information? 

Both above questions have been addressed in the momentum fraction representation and, more 
recently, in the framework of the 'dual' GPD parameterization SU ED]- We shall recall the 
known partial answers to those questions in the momentum fraction representation and add a 
rarely recognized GPD aspect, namely, duality [121 [31], with the focus on the phenomenological 
applicability. Before we give the technicalities in the main body of the paper, we briefly outline 
the answers here. In our opinion it does not matter in which representation they are formulated, 
and in this paper we prefer the momentum fraction representation. 

The answer to the first question arises from analyticity and it has been clearly spelled out for 
a fixed resolution scale in Ref. [13]. Let us recall that the DVCS amplitude is considered to be a 
holomorphic function, see, e.g., Refs. [1U SSI SSI EH SHI HZ] , whose real part can be expressed in 
terms of its imaginary part, i.e., its s- and w-channel discontinuities, via a single variable dispersion 
relation for fixed t. Moreover, to LO the imaginary part is given by the GPD on its cross-over 
trajectory (rj = x), which separates the meson- like t-channel exchange interpretation of the central 
region (\x\ <r]) from the partonic s-channel exchange view of the outer region (rj < x < 1). Hence, 
apart from the cross-over trajectory the GPD for fixed resolution scale cannot be experimentally 
accessed. The dispersion relation gives us a considerable handle on the GPD at the cross-over 
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trajectory and constrains its functional form, even if the modulus and phase of the amplitude are 
measured only in a restricted phase space. Therefore, the answer to the first question guides us to 
a simple phenomenological concept in which the dispersion relation is used to pin down the value 
of the GPD on the cross-over trajectory. 

Suppose we know now the GPD at its cross-over trajectory from an (ideal) measurement for 
fixed photon virtuality Q 2 , as illustrated below in Fig. [U The deconvolution problem, which then 
pops up, is how the GPD at 77 = x is related to its value at any 77, i.e., for the whole support. 
Lorentz covariance guarantees that a GPD possesses a duality property [121 [M], i.e., it can be 
mapped from the outer region into its central one, where the GPD is continuous on the cross- 
over trajectory. Obviously, one can draw a priory any surface over the outer region of the GPD 
support that connects the curves given by the GPD at 77 = x and 77 = for < x < 1, see Fig. [TJ 
Positivity constraints, see Ref . [4*8l |4*9] and references therein, lead to constraints on these surfaces. 
The dual surfaces in the central region are then governed by Lorentz covariance. Each of these 
surfaces corresponds to a GPD model, and so the deconvolution problem has infinite solutions. 
The problem with present ad hoc GPD models is that usually one solution is picked out, with 
parameterization of the GPD at the cross-over trajectory being too rigid. We stress that evolution 
of GPD with change of the resolution scale provides a supplementary handle on the deconvolution 
problem, as advocated in Ref. [51] . However, for a small lever arm in the photon virtuality, 
only if some holographic principle exists, rigidly relating shape of this surface to values along the 
77 = x cross-over trajectory, the deconvolution problem can be strictly solved. But even if this 
hypothetical principle remains unknown, one still has various theoretical and phenomenological 
tools at hand to address the deconvolution problem and, in particular, the related problem of 
decomposition of the quark angular momentum. 

In this paper we introduce a new set of such tools, namely GPD sum rules, and we illustrate 
how an alternative phenomenological approach to DVCS might look like. The outline is as follows. 
To shed light on unrealistic features of the GPD models, presently used in phenomenology, we give 
in Sect. [2] a short overview of experimental findings and their confrontation with GPD models, 
mainly for the DVCS process. The reader who is not interested in this phenomenological aspect 
can safely skip this section. In Sect. [3] we employ previous results [121 [M] to clearly spell out that 
GPDs have an internal duality, relating the outer and central regions. We derive then various 
GPD sum rules. In Sect. H]we discuss the phenomenological applications of these sum rules and 
also give a preliminary example for fixed target kinematics. Finally, in Sect. [5] we summarize our 
findings and conjecture that GPDs should be governed by a holographic principle, which might 
be naturally related both to the duality interpretation of strong interaction phenomena, reviewed, 
e.g., in Ref. [52], and to the AdS/CFT conjecture PI [MIES]. 
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2 Mainstream GPD models versus experiment 



Numerous experimental, theoretical, and phenomenological efforts, reviewed in Refs. 0, [7J , have 
been spent during the first decade of GPD phenomenology in an attempt to go beyond the basic 
GPD ansatz, given by the product of PDF and elastic form factor. Fortunately, there is tremendous 
progress on the experimental side and, certainly, some theoretical understanding is reached, too. 
Unfortunately, some rather hypothetical elements became an important part of the mainstream 
phenomenology, while at the same time some well-understood theoretical features are incompletely 
implemented. In what follows, we would like to make it clear that the observed phenomenological 
problem, namely, that certain GPD models fail to describe the DVCS data [21 [23 ESI ED], indicates 
only that the ad hoc GPD model approach in its present form is an inefficient tool for the future 
and that any general conclusion about the GPD formalism drawn from the failure of certain GPD 
models should be considered as a speculation. Below we briefly describe the situation and attempt 
to separate different problems, which are often mixed up. 



2.1 GPD basics and terminology 

GPD models, used in phenomenology so far, are mainly based either on a spectral representation 
for the Green functions of light-ray operators, so-called double distributions (DDs) [U [57], or on 
the collinear conformal partial wave expansion^ [59] , known in various versions [60l |6T| |62| 1631 
H2J EH E5]- Since they all by definition represent the same field theoretical object, any GPD 
model, as long as it respects basic properties, such as Lorentz covariance, can be represented in 
any of the named representations. Another popular GPD representation is given as overlap of 
light-cone wave functions [661 EH EH] . We emphasize that the one-to-one correspondence of the 
DD and overlap light-cone wave function representation [611 E3 EE] should hold on general field 
theoretical grounds, too. This has been demonstrated within wave function models in Refs. [691 
[701 ITT] . In the overlap light-cone wave function representation the positivity constraints [HI H9] 
are explicitly satisfied^] These constraints are exactly valid at LO approximation; however, their 
power drastically diminishes if one calculates the physical amplitude. 

2.1.1 Double distribution representation 

The DD representation arises directly from diagrammatical considerations [HE?]. We consider it 
to be the most general representation that inherits properties of field theory. Here, the partonic 



lr The group theoretical aspects and QCD applications of the collinear conformal group are given in Ref. [58] . 
2 The solution of these general positivity constraints can be interpreted as a two body light-cone wave function 
overlap representation for a collective spectator |50j . 
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quantum numbers are the momentum fractions of two partons. We recall that the first sophisti- 
cated GPD ansatz was elaborated in Refs. [721 173] for t = via a DD ansatz, which factorizes in 
a PDF and a single variable profile function whose shape, viewed as a meson distribution ampli- 
tude, was assumed to be convex. For brevity, we call this the Radyushkin's DD ansatz (RDDA). 
Because of its convenient parameterization, it was adopted by many model builders in various 
t-decorated versions and within generalizations of the profile function; however, the convexity of 
its shape was not touched upon. In particular, the RDDA [72], [73] has been utilized in Ref. [7J] to 
build up a GPD model, guided by the results of the chiral quark soliton model (xQSM), decorated 
with a factorized or Regge inspired t-dependence, and based on certain PDF parameterizations. 
This GPD model is now numerically implemented in the popular VGG code Here a descen- 
dant of the D-term has been utilized to parameterize the GPD E in terms of the quark angular 
momentum. We follow the common terminology and call this code the VGG model. We note 
that VGG results might be different from other GPD models that were written down in a similar 
fashion within RDDA, e.g., in Ref. [2"7j . 

It was already emphasized in Refs. [721 [73] [75] that more realistic GPD models should possess an 
intricate interplay of momentum fraction and momentum transfer squared dependence. It turned 
out that guidance for a more realistic t-dependence comes from field theory inspired models. 
Such an improvement is now implemented in GPD models for the pion [BH] arid nucleon [70] [7T] . 
All these investigations started from overlap representations with power-law wave functions, and, 
aiming at different goals, they finally provided rather analogous expressions for the resulting DDsJj 
Such models can also be obtained from the evaluation of (general) triangle Feynman diagrams. 
To emphasize the physical content of these GPD models, one might also denote them 'collective 
spectator quark models' Q For both t] = and r\ = x the t-dependence of these models dies out in 
the limit x 1; see, e.g., discussion in Ref. [71]. The behavior in the former case coincides with 
lattice measurements of generalized form factors [TTJ [78] [79] [80] [81] [82] [83]. A more general DD 
ansatz in the form of a factorized RDDA with the common improvement of Regge behavior has 
been written down for certain proton GPDs in Ref. [71]; for the pion case see Ref. [HH]- We have 
also convinced ourselves that one should not insist on a convex shape of the profile function, and 
could, e.g., add 'higher' modes, to get rather flexible GPD models. In this way one can replace 

3 It was observed in Ref. [69] that the functional form of light-cone wave functions cannot be arbitrarily chosen 

to derive the DD representation. Therefore, the authors of Ref. [70] started with a covariant quark model to obtain 

GPD models, formulated in DD representation. In Ref. [7T] it was stated that in the light-cone wave function 

overlap representation the longitudinal and transversal degrees of freedom are tied to each other due to Lorentz 

symmetry. Hence, arbitrary modelling of light-cone wave functions will usually violate the polynomiality property 

of GPD moments, see, e.g., Ref. [76] . 

4 There is a huge literature in which GPDs have been studied within dynamical model approaches. The reader 

may find references in Refs. [61 [7T1 [7S] , 
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ad hoc GPD model approach with a more flexible fitting procedure. 
2.1.2 Conformal partial wave expansions 

The expansion of the GPD in collinear SL(2, M.) conformal partial waves arises naturally from 
the solution of the evolution equation to LO accuracy. The parameterizations, applied in phe- 
nomenology, are named 'dual' [M] or Mellin-Barnes representation [12]. In our understanding, 
the name 'dual' primarily refers to the use of crossing symmetry [H (S3 ES|, which is, by itself, in- 
dependent of the chosen representation. For the sake of clarity, we follow the authors of Ref. [M] 
and denote their representation 'dual' in quotation marks. The term dual (without quotation 
marks) we use for a general GPD property, as used above and explained further in Sect. [3j This 
concept was partially utilized for evolution kernels [S3 d] and has been adopted to GPDs [32"] . 
supplemented by a constructive proof [34j . As emphasized in Ref. [12] it is also tied to crossing, 
again independently of the representation. The term Mellin-Barnes representation refers to a 
specific integral representation of the collinear conformal partial wave expansion, closely related 
to the well-known Sommerfeld- Watson transform, in which internal GPD duality is manifestly 
incorporated, cf . Ref. [HU ES] • 

In the Mellin-Barnes representation, besides the complex-valued conformal spin, the angular 
momentum of the t-channel SO (3) partial wave might be employed as a label. In the 'dual' 
parameterization a forward-like momentum fraction (conjugated variable to the conformal spin) 
and the difference of integral conformal spin and t-channel angular momentum appear as labels. 
This makes contact to hadronic physics [H], in particular to Regge phenomenology [31]. In the 
minimalist version of this model the expansion of conformal moments in terms of SO (3) t-channel 
partial waves is restricted to the leading one. The inclusion of the next-leading SO (3) partial wave 
has been called minimal 'dual' parameterization [HI]. More recently it has become clear how the 
inclusion of non-leading SO(3) partial waves can be effectively achieved in the 'dual' parameteri- 
zation |41| . The minimalist and minimal version in the 'dual' parameterization corresponds to the 
leading and next-leading SO (3) partial wave approximation. We use the terminology: leading (i.e., 
minimalist) SO(3) partial wave and next-leading (i.e., minimal) SO(3) partial wave approximation 
|34j . GPD modelling in the 'dual' and the Mellin-Barnes integral parameterizations is presently 
done with different phenomenological emphases. In the former parameterization the model aspect 
with respect to the approximation of SO(3) partial waves, see Sect. I4.1[ is pronounced [101 HI], 
while in the latter, up to now applied only to small XBj physics, the Regge aspect [31] and the 
resummation of SO (3) partial waves [88] is given priority. 
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Figure 1: The left panel displays the spectator model of Ref. [7T], as specified there in Eq. (25), 
with t = t m in — 0.25 GeV 2 . The red solid curve visualizes the experimentally accessible information 
for fixed photon virtuality, i.e., the GPD at the cross-over trajectory rj = x, while the dashed one 
is an example curve that can be used to evaluate the real part of CFFs via Eq. (j3J). In the right 
panel the same GPD at r\ = x (solid) is displayed versus the Bjorken-like variable X = 2x/(l + x) 
and it is compared with a VGG-like (dashed) and the leading SO (3) partial wave approximation 
of the Mellin-Barnes (dash-dotted) ansatz. Both ansatze are based on the zero skewness spectator 
GPD model and possess a factorized t-dependence. 

2.1.3 Other GPD expansions 

One might also introduce an expansion with respect to a complete basis of orthogonal polynomials, 
such as Gegenbauer [59] or Bernstein [89] Such expansions can in principle be expressed in the 
two considered ones. The expansion in terms of Bernstein polynomials has been proposed once in 
deep-inelastic scattering (DIS) [90]; however, it could not compete with the standard approaches, 
employed in global DIS fits. Also, we know from experience that the expansion of GPDs with 
respect to a complete basis of orthogonal polynomials turns out to be numerically inefficient at 
smaller values of x-q] [92] . 

2.2 Phenomenological applications 

The phenomenological application of GPDs for the deep virtual leptoproduction of mesons or 
photon is based on factorization theorems, proven to leading order in the 1/ Q expansion [4], 15) . 

interestingly, in Ref. |91j the authors started from a covariantly formulated spectator model; however, consider 
Lorentz covariance (polynomiality of GPD moments) as a DD hypothesis. Finally, they are doing GPD modelling 
for rj = 0. To include skewness dependence, they parameterize the GPD in terms of a few Bernstein polynomials 
and rely on lattice data [53] ■ 
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Thereby it is assumed that GPDs are continuous at the cross-over point. It is commonly assumed 
that the GPD framework is quantitatively applicable for DVCS at a scale of the order of one or few 
GeV 2 , while for meson production it might be of the order of 10 GeV 2 . Unfortunately, only very 
little is known about the size of power suppressed contributions and it remains an open problem 
to pin down the onset scale for the GPD formalism. This problem might also be addressed from 
the phenomenological side. In the following we assume that the GPD formalism is quantitatively 
applicable to leading order of the perturbative expansion in the QCD coupling constant for DVCS 
measurements in JLAB, HERMES and Hl/ZEUS kinematics, where the photon virtuality is larger 
than 1 GeV 2 . 

The specific GPD models are all more or less constrained by the same phenomenological input, 
PDFs, form factors, and sometimes model or lattice estimates. Roughly spoken, the difference 
in the physical content of various GPD models at t] = x is more or less accidental and arises 
from simplifications, mainly done for convenience, in a given representation. Obviously, the model 
dependence is still visible in the described deep virtual amplitudes; however, the parameterization 
thereof is effectively not under control. This is illustrated in the right panel of Fig. (TJ where we 
show the GPD at its cross-over trajectory within the spectator model (solid) of Ref. [71], the 
leading SO (3) partial wave approximation in the Mellin-Barnes representation (dash-dotted), and 
a VGG-like ansatz (dashed). For the latter two models we have used a factorized t-dependence, 
and the form factor and PDF were evaluated in the spectator model. There is one rule which 
GPD model builders like to follow; namely, they mostly avoid invisible degrees of freedom (those 
that die out in the forward limit). In other words, one relies on the implicit assumption that the 
holographic principle holds true. Exceptions to this rule are the pion pole contribution [9"3"II71] and 
the so-called D-term [HI], where the latter is entirely related to the subtraction constant in the 
dispersion relation [33]. In the SO (3) t-channel partial wave expansion both of these exceptions 
are assigned to J = contributions, as shown for the subtraction constant in Sect. 131 In the 
context of the parton model jHSj, and within its field theoretical model realization for forward 
Compton scattering [96J, such a J = term relates low- and high-energy physics. We emphasize 
that there are various possibilities to embody these J = terms into GPD models, see, e.g., 
Refs. [931 E7J ESI [701 [331 EE], and they might lead to rather different ad hoc GPD 'predictions'. 

2.2.1 Ad hoc GPD phenomenology 

The qualitative failure of the VGG model in describing virtual electroproduction of p° mesons in 
the resonance, i.e., large xb], region [25] is partly related to its ad hoc t-dependence. Certainly, here 
one has to be concerned about the applicability of perturbative QCD (and, by the way, also about 
the onset of the Regge regime) and so one might expect in this region both large radiative and 
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power suppressed corrections. However, since the Feynman mechanism is somehow incorporated 
in the GPD formalism, we might expect that the t-dependence of the hadronic amplitude is 
qualitatively embodied in the GPD; see Ref. [981 199] for a perhaps more appropriate realization 
of the GPD formalism in this kinematics. The t-dependence of the cross section diminishes in 
the large xb] region [25], which is consistent with the aforementioned spectator GPD model |71jr 
The notorious normalization problem of the cross section in this kinematical region is to LO 
approximation entirely related to the end-point behavior of the VGG model. It was cured in 
Ref. [25] by a D-term inspired addendum, concentrated in the central GPD region and vanishing 
in the forward kinematics. We will also illustrate in Sect. 14.21 that one can easily construct GPD 
ansatze which could fix the normalization problem in a different manner. In fact, it is demonstrated 
in Fig. [U(b), that a simple spectator model (solid curve) possesses in comparison to a VGG-like 
ansatz (dashed curve) an enhancement effect in the large momentum fraction region. 

The DVCS amplitude is accessible in the deep virtual electroproduction of photons, where due 
to the interference of the bremsstrahlung and DVCS processes one can also reveal the phase of the 
DVCS amplitude. The VGG modej^ also fails to describe the beam spin asymmetry and polarized 
cross section for unpolarized proton in the fixed target JLab kinematics, see e.g., Ref. [13], since 
here it does not offer the possibility to decrease the normalization of the resulting amplitudes, 
compared to the basic GPD ansatz. 

The normalization problem was already studied in Refs. [23 EE] and to some extent cured 
by reducing the amount of sea quark contributions (introducing so-called model A or B, which, 
however, have an unrealistic t-dependence) . In fact, adjusting the normalization of the DVCS 
amplitude by suppression of sea quark contribution^ makes it possible to describe existing DVCS 
data, related to beam spin asymmetries or polarized cross section measurements. Even semi- 
quantitative discussions [271 12H] seem to provide empirical understanding, when employed within 



6 We thank the authors of Ref. asking the question "GPDs or not GPDs?", for clarification. Their conclu- 
sions about GPD models within DD representation refer only to specific models, considered there. In particular 
the statement "the very flat t slopes observed can only arise from GPD contributions not constrained by the FF 
[form factor] sum rule" is not meant to be valid in general. Indeed, simple spectator GPD models, formulated in 
DD representation and satisfying the form factor sum rule, explain the observation. 

7 We would like to emphasize that the 'easy to tune' options in GPD codes indirectly suggest the possibility 
of easy phenomenological GPD interpretation; however, we fear that this can be misleading. For instance, it was 
found in Ref. [20] that the D-term within the VGG model is disfavored from the HERMES DVCS measurements. 
Since the xQSM [74| . 1100] and now also lattice calculations, e.g., in Ref. [83], indicate that the D-term is sizeable, 
our interpretation of this statement is that the VGG model, the incorporation of the D-term equivalent part (see 
Sect. 3.3) in the model or a combination of both, is improper. Also, the model dependent constraints on the quark 

angular momentum in Ref. [20] illustrate that the model dependence is not under control. 

8 Sea quark contributions and in our understanding also valence singlet quark ones are neglected in the DVCS 
results of Ref. . 
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the impulse approximation in a prediction of beam spin asymmetry for DVCS on nuclei [101j . If 
this 'prediction' turns out to be in disagreement with increasingly precise experimental data, it 
will be exciting to employ the GPD formalism as a tool to explore the interplay of nucleon and 
partonic degrees of freedom beyond the impulse approximation. 

The minimal(ist) version of the 'dual' parameterization [M] has been confronted with DVCS 
data in Refs. [29] and [30] and describes asymmetries and polarized cross sections. Finally, all 
these GPD models fail to describe the unpolarized DVCS cross section, measured by the Hall A 
collaboration at JLab |14j . 

The understanding of GPDs in collider kinematics of the HI and ZEUS collaborations [9j 
[T2"| [TBI [T5] was in poor state, too. Claims that the DVCS cross section can be described within 
published GPD models to LO accuracy are based on either a problematic implementation of the t- 



dependence [271 [29] PJ violation of the perturbatively predicted scale change [29[|102j J 10 l or violation 
of basic GPD properties [T02]F 

We stress that the physical content of the LO normalization problem was known to the authors 
of Ref. [102] . Although they started to analyze the model problems within the RDDA, they finally 
gave up on the DD representation. Armed with the aligned jet model, a reasonable physical 
picture, however, violating Lorentz covariance, they came up with an inconsistent GPD model 
(see footnote [TT]) , which must be rejected. After restoring Lorentz covariance in their model one 
simply finds a PDF which is decorated with t-dependence — a special case of the RDDA. 

The normalization problem can be naturally overcome by inclusion of radiative corrections 
as advocated in Refs. [311 [32]. However, it was at the same time observed that GPD ansatze 
within different PDF parameterizations, describing the same set of unpolarized DIS data, result 
in rather different DVCS predictions [92] • This is quite understandable if one remembers that 
already the description of the unpolarized DIS structure functions is a fine tuning problem. There 
the perturbative expansion of the evolution is unstable in the small xej region, see discussion 
in Ref. [103] , and the same thing now happens in the DVCS case [31]. We emphasize that 
this does not disqualify the perturbative approach to DIS and DVCS at small x^, because the 
evolution operator is universal, i.e., process independent. Also a resummation procedure might 



9 Since the ^dependence of the sea quark GPD is unrealistic, we rejected the model of Ref. [27] for some time. 
However, it describes various DVCS data and we employ it for error estimates in Sect. 4.2. The description of data 

in Ref. [29] arises from a GPD model that induces an unrealistic shrinkage effect of the diffractive forward peak. 

10 The incorporation of the Q 2 -dependent i-slope (32) in Ref. [102] . also adopted in Ref. [3S], violates the LO 

scale change prediction. 

11 The aligned jet GPD model [102J is based on a truncated Taylor expansion in rj for all GPD moments. The 
inconsistency of the proposal [102] was revealed in Ref. [46]. According to Eq. (I2ip . the moments in the 'full' 
model, evaluated from Eq. (21, [102] ). are ry-independent and the resulting GPD is the PDF decorated with a 
^-dependence. 
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Figure 2: The left panel shows a preliminary LO fit result [88J to the HERA DVCS data set versus 
the photon virtuality Q 2 : ZEUS p] with W = 89 GeV (triangle, solid curve) and HI with W = 
82 GeV within the HERA I pj] (empty circle, dash-dotted curve) and HERA II [18] (full circle, 
dash-dotted curve) run. In the right panel we show the skewness ratio R, defined in Ref. [18] as 
observable, (dashed curve) to LO approximation, i.e., R ~ H(x,x,t = 0, Q 2 )/H(2x, 0, t = 0, Q 2 ) 
for fixed W = 82 GeV, and the ratio r = H{x, x,t = 0, Q 2 )/H{x, 0,t = 0, Q 2 ) for fixed x = 1(T 3 
(solid curve). Note that R « 2"(°> s2 ) r with a(0, Q 2 ~ 4 GeV 2 ) ~ 1.2. 

stabilize the perturbative results; for a discussion see, e.g., Ref. [104] . In particular, it will be 
phenomenologically important for a perturbative GPD analysis of the deep virtual vector meson 
electroproduction, see Ref. [1U5] for details. 

2.2.2 A fitting procedure 

The ad hoc GPD ansatz approach is surely not a proper framework for revealing the partonic 
content of the proton. Insight in the perturbative GPD approach is needed, too. We started 
to explore an alternative framework for DVCS in the small kinematics |106[ [341 [88] up to 
next-to-next-to-leading order |1U7[ I1U81 11U9[ II 101 [53] and would like here to present what we have 
learned. Since radiative corrections to the DVCS amplitude (or DIS structure function Ft) can 
be essentially absorbed by a scheme change [27], the description of experimental DVCS data at 
LO accuracy is the indispensable key to a GPD understanding. 

In the small x-Qj kinematics the functional form of the amplitude simplifies drastically and the 
reggeized t-channel view gives guidance for modelling. This is a pure empirical recipe, where one 
has to replace the classical pomeron trajectory by an effective one. The effective pomeron intercept 
is an integral part of the universal GPD and it is considered to be skewness independent, i.e., the 
GPDs or DVCS amplitude have the same intercept as the corresponding PDFs or DIS structure 
function, respectively. This is supported by both a diagrammatical t-channel ladder analysis [33] 
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and the experimental data. 

For a fitting routine it is crucial to have full control over the normalization of the DVCS 
amplitude. This can be achieved within invisible terms that die out in forward kinematics. The use 
of conformal symmetry enabled us to study the GPD approach up to next-to-next-to-leading order 
in perturbation theory |110ll33| IM]. The numerical advantages of the Mellin-Barnes representation 
then allow for a simultaneous fitting procedure of DVCS and DIS data. The essential lesson from 
our GPD fits is that a leading SO(3) partial wave approximation (minimalist 'dual' model) does not 
work in the small xb] regime at LO, but only at perturbative next-to-leading order or beyond [M] . 
Within the conformal partial wave expansion, a more flexible parameterization of GPDs has been 
worked out by several researchers; however, it is only partly published |111U1121HI)II8"8"] . After such 
an improvement, the experimental data for collider kinematics can be fitted at LO accuracy either 
within the inclusion of the next-leading SO (3) partial wave [11 Xj or a flexible model dependent 
SO(3) partial wave resummation [88]. A preliminary fitting of results to the DVCS cross section, 
measured by the HI [T3"l [T8"] and ZEUS [12] collaborations, versus the photon virtuality Q 2 [SB] is 
displayed in the left panel of Fig. [21 Our result is that the dominant sea quark GPD at low x shows 
to LO accuracy almost no skewness effect [88]; see the solid curve in the right panel of Fig. [2j We 
should here give attention to a different GPD interpretation of the same DVCS data. Namely, 
the measured ratio R [IB] and our corresponding GPD fit (dashed line), might be viewed as a 
large skewness effect. Actually, using a better suited r function, we reinterpret this as (almost) 
no-skewness effect. Such a vanishing skewness effect was for a lower resolution scale advocated 
in Ref . [H] , but it turns out that it holds for the experimental lever arm of ~ 3 ■ • • 80 GeV 2 
|18j . In the published GPD models this no-skewness property was either never present or, if it 
was taken at low input scale [H], the ansatz spoils it due to the evolution [H]. Unfortunately, 
the large evolution effect was propagated as a 'prediction' [113] , it is 'naturally' implemented in 
approximated versions of conformal partial wave GPD representations, and it occurs in a weaker 
form for RDDA, see Ref. |114j . We add that in the phenomenology of diffractive vector meson 
production it is often taken for granted that the GPDs possess a positive non-zero skewness effect, 
e.g., in Ref. [115] . 

2.2.3 Lessons from failure and success: a speculative GPD picture 

From a phenomenological failure of present GPD models one can mainly conclude that its ad hoc 
realization was improper. It is senseless to criticize GPD representations in general, which will 
just lead to a step backwards. GPD features which have been understood in one representation 
can be incorporated in the other representation, too. Furthermore, from the confrontation of 
experimental findings and present GPD models one cannot provide a reliable partonic interpreta- 
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tion, in particular for the D-term equivalent part and the quark angular momentum. The critique 
of the perturbative approach offered so far from the outcome of such a confrontation should be 
considered as speculation. In particular, from what was said above for DVCS, it follows that a 
perturbative GPD analysis for deeply virtual vector meson electroproduction at small x-Qj also has 
to be reconsidered and extended to next-to-leading order. 

The failure of the mainstream ad hoc GPD models should be an incentive to learn the lessons 
and find improved versions. This has been achieved to some extent with respect to the interplay 
of momentum fraction and t-dependence [69j [ZQl EI] - Pragmatical solutions to fit data within the 
GPD formalism are obvious, namely, one has to exploit degrees of freedom that are invisible in the 
zero skewness case, e.g., more flexible profile functions. This sounds simple; however, it requires 
some technical effort and, more importantly, GPDs are intricate functions and so one can easily 
get 'lost'. Even in the in some sense much simpler case of small XBj kinematics, three different 
groups of authors attempted in a period of five years or so to go for a GPD interpretation of the 
ZEUS/HI DVCS data, where two of them gave up basic (and realistic) GPD properties. In our 
opinion the message from that for the analysis of fixed target experiments is that one should first 
utilize dispersion relation techniques, to pin down the GPD on its cross-over trajectory. 

Nevertheless, we attempt in the following to catch the possible qualitative features of a 'more 
realistic' GPD at its cross-over trajectory from the phenomenological findings and we would like 
to share our speculations on how it might look, away from the cross-over trajectory. These un- 
orthodox considerations are disputable and we consider the outcome only as a working hypothesis. 

Let us briefly summarize the phenomenological findings. Present ad hoc GPD models, perhaps 
apart from the VGG one, are in the valence region consistent with (beam) spin asymmetry and 
polarized cross section measurements. If they are tuned to experimental data they might possess 
some kind of predictive power for certain observables. There is now a clear understanding of the 
dominant sea quark GPD in the small xbj region, namely, the skewness effect to LO accuracy is 
compatible with zero. The extension of such an ansatz to larger values of x induces a wrong scale 
dependence. The large size of the unpolarized DVCS cross section in the valence region is hardly 
understandable within present GPD models. The virtual p° electroproduction data [25] in the large 
XBj region are consistent with respect to the t-dependence of simple spectator models. A precise 
quantitative understanding of the cross section is for various reasons rather challenging. We add 
that the employment of unitarity constraints in a Regge inspired model analysis [SB] indicates that 
for DVCS already the valence region cannot be correctly described with a few reggeized t-channel 
exchanges. 

Since the real part of the DVCS amplitude can be determined from a dispersion relation, 
a GPD ansatz at the cross-over trajectory must simultaneously give the correct description for 
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the imaginary part of the amplitude in the whole phase space. This is certainly not the case 
for present GPD models; however, this GPD aspect can be straightforwardly analyzed using the 
growing amount of DVCS data. 

On qualitative grounds we conclude from the experimental findings that a GPD on its cross- 
over trajectory is at small x governed by the reggeized t-channel view, where the skewness effect 
might be zero, while at large x there is a large skewness effect and the t-dependence diminishes. 
Such a behavior is supported by field theory inspired spectator models (69J E2 d]; see, e.g., the 
discussion in Ref. [71], where, however, the Regge behavior at small x is absent. Implementing 
Regge-behavior from the s-channel view, along the lines of Ref. [116j . yields at t = a Regge 
'improved' RDDA that possesses an experimentally unfavored large skewness effect at small x. 
To find a more realistic behavior it is perhaps more appropriate to have the t-channel view by 
considering the diagrammatical ladder, as started in Ref. [H] . That this is a reliable concept has 
been demonstrated for PDFs, see, e.g., Ref. [117] . 

We finally mention that due to the evolution the GPD value in the outer region is reduced 
and 'flows' into the central region to 'fill' it, i.e., yielding a growing value of the GPD. From the 
experimental observation that the skewness effect at very small x 'generically' does not depend 
on the resolution scale, one might draw a speculative conclusion. Namely, loosely spoken, the 
'flow' into the central region has already stopped and the GPD in this region is mostly 'filled' at a 
resolution scale of about 2 GeV or perhaps even less. This seems to fit the qualitative experimental 
findings at large a?Bj, indicating a 'big' GPD contribution in the central region [25]. Utilizing 
the evolution equation, this pictorial speculation can be exactly formulated as a mathematical 
problem; however, this is beyond the scope of the paper. 

3 GPD duality and families of GPD sum rules 

In this section we deliver the tools that allow for a straightforward GPD analysis of the photon 
electroproduction data in fixed target kinematics. Thereby, we adopt well-known techniques, 
utilized in the phenomenology of on-shell processes and combine them with the duality property 
of GPDs gang]. 

3.1 Preliminaries 

To address the problems, we consider virtual Compton scattering off a proton in the generalized 
Bjorken region [I]: 

7*(<7i) P(Pi) -> T W (?2) p(p 2 ) , Q 2 = + <? 2 ) 2 /4 -> oo , (1) 
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where the momentum transfer squared t = A 2 = (p 2 — Pi) 2 and the two scaling variables 

t Q 2 A-g 
P ■ q P ■ q 

are fixed. Here P = p\ + p 2 , A = p 2 — pi and q = (q\ + <?2)/2. We shall always take the values 
of £ and 77 to be positive. In the twist-two approximation the crossing- invariant (i.e., for s *-* t) 
photon asymmetry parameter is given by their ratio: 

The amplitude is parameterized in terms of the Compton form factors (CFFs), which in pertur- 
bative LO, Born or hand-bag approximation, read for leading twist-two: 



0, t, Q 2 ) L =? / dx ( - 1 - =f - ; 1 , ) F(x, r, = t, /i 2 = Q 2 ) . (4) 



X — x — ie £ + x — ie 

Here the scale fi 2 is a factorization scale at which the hard-scattering part and GPDs are factorized 
and it is equated to the characteristic scale Q 2 of the process (the common choice). The explicit 
dependence of the hard-scattering part on /i 2 appears at next-to-leading order. In our shorthand 
notation CFFs will be denoted T and the GPDs F, cf. Ref. 



^={~'|}(^,Q 2 ), F=||'|}(x,r / = ^,t,/, 2 ) forparity j^} . (5) 

The sign =)= on the r.h.s. of Eq. (j3j) refers to the symmetry with respect to x. This convention will 
be consistently used below without an explicit reference to it. Note that antisymmetry (symmetry) 
with respect to x corresponds to parity even (odd) CFFs, which also have signature +1(— 1). This 
definite signature originates from symmetry under the exchange of s- and it-channel. Here we did 
not indicate fractional charge squared factors, which can easily be restored by replacing F with 
Tlia=u d ■■■ Qa a F- The GPDs F(x, rj,t, Q 2 ) are real- valued functions that are symmetric in r\ [118] . 
and their jth Mellin moments are polynomials in 77 of order j + 1 (polynomiality property). To 
be on safer ground concerning the analyticity hypothesis, we consider in the following only the 
Euclidean region, i.e., < 1. This is ensured by the condition £ > rj. Unfortunately, apart from 
the contribution to virtual two-photon exchange corrections |119j . the CFFs are only accessible 
in DIS (0 = 0, i.e., rj = 0) and DVCS (0 = 1, i.e., rj = £). In the former case we extract PDFs, 
equal to GPDs with 77 = and t = 0, and in the latter case we find the GPDs on the cross-over 
trajectory x = rj. For pedagogical reasons, however, we shall first do a gedanken experiment and 
suppose that the Compton amplitude can be 'measured' in the whole Euclidean region. 

The obvious part of the answer to the question 'What GPD information can be revealed from 
experiment?' is that the imaginary part of the CFF (j3J) is given by the GPD combination: 

FT( X , 77, t, Q 2 ) = F{x, v, t, Q 2 ) =F F(-x, 77, t, Q 2 ) = -%nLF{t = x, = rj/x, t, Q 2 ) . (6) 

7T 

1(3 



This formula tells us that by varying the photon asymmetry parameter, i.e., < $ < 1, the GPD 
is scanned in the outer region 

0<77<2C = ^<1. 

This is analogous to DIS, where the structure functions are to LO equal to PDFs. The r.h.s. in 
Eq. ([6]) and so the GPD F T (x,r],t,Q 2 ) with rj < x might be viewed as the spectral function 
for this Compton scattering process, given by the s-channel cut. Hence, the outer region of the 
GPD F(x, t], t, Q 2 ) with 77 < x is expressed by the difference or sum of s- and w-channel cuts. In 
accordance with the standard definition of the PDFs, the GPD is understood for 77 < x as the 
combination F T (x > rj, rj, ■ • • ) = [ q F + ^F] (x, rj, • • • ) of quarks and antiquarks, where 

\«F) K ' 7 \=fF(-x,77, •••)/ \antiquark ~' v; 

The definite symmetry of F T GPDs under reflection x —* —x allows us to extend them uniquely 
to negative values x < —77. The support extension into the central region x G [—77, 77] should be 
done in such a way that the polynomiality condition is satisfied. As we shall see in Sect. !3.3l below. 
analyticity guarantees that this procedure is unique. 



3.2 Sum rules 

To get deeper into the problem, we consider the CFFs as holomorphic functions of Mandelstam 
variables which satisfy the Schwartz reflection principle. Furthermore, we assume that for fixed 
(negative) value of t in the Euclidean region the CFFs possess only elastic poles, as well as s- 
and u-channel discontinuities on the real axis. Hence, CFFs satisfy a single variable dispersion 
relation, where t and the photon virtualities are considered as fixed variables; see, e.g., Refs. (HI 
SSI H31 [3U H6] . Taking into account that they have a definite signature, this can be equivalently 
rewritten in terms of the crossing-symmetric variable v = (s — u)/4M, related to the energy of 
the initial photon, or in terms of the scaling variable £ = Q 2 /2Mu, originally defined in Eq. (j2J): 

:F(£, 0, t, Q 2 ) = - [ d£' ( I , =f I . ) 3mf(e' - iO, 0, t, Q 2 ) + C^, t, Q 2 ) . (8) 
^ Jo £ + £' -lej 

For convenience, we here set the upper integration limit £ cut = £ po i = 1/ (^ + t/AQ 2 ) to 1, according 
to the Bjorken limit, so that the elastic pole is included in the imaginary part, rather than in the 
subtraction constant. We emphasize that the correct upper bound can easily be restored by the 
substitution 1 — > 1/(1 + t/AQ 2 ). The subtraction constant is fixed by the value of the CFFs at 
£ = 00. It appears in the CFFs 7i and S, vanishes for the CFFs 7i and S, and is perturbatively 
predicted to be zero for the combination 7i + S, see Ref. [31]. Thus, we can define: 

C(0, t, Q 2 ) = C e (0, t, Q 2 ) = £(£ = 00, t, Q 2 ) , C n (#, t, Q 2 ) = -C{0, t, Q 2 ) , (9) 



17 



where = Cg = 0. 

The dispersion relation (jSJ) states that both the imaginary and real part of CFFs contain, up 
to a possible subtraction constant, the same information: 

&e.F(f , 0, t, Q 2 ) = ^PV J d£ L-^j, =F — ^) - zO, 0, t, Q 2 ) + Cjr(0, t, Q 2 ) . (10) 

Let us loosely have a look at this statement without discussing any mathematical refinements. 
After extending the support of SmjF(^, . . .) to |£| > 1 and using the symmetry property, we 
express Eq. (I10p as a (singular) Fredholm type integral equation of the first kind with constant 
limits (±oo). Briefly speaking, the real part of . . .) — C(. . .) is the Hilbert transform of its 
imaginary part. The inversion reads then 

i r°° l 

sm^(e, 0, *, g 2 ) = -pv / dC —— [&eF(t', 0, t, g 2 ) - cM0, t, g 2 )] . (11) 

Alternatively, by again making use of the spectral condition QmjF(£, . . .) = for |£| > 1, we can 
take the solution of the Fredholm type integral equation with constant limits of integration (±1); 
e.g., from Ref. [T20] : 

^(e,^,t,g 2 ) = ^Hgpv r^^ ^ene^g 2 ) for (12) 

In conclusion, assuming the specific analytic properties, we can state that the imaginary and 
real part of CFFs, having definite signature and a vanishing subtraction constant, have a one-to- 
one correspondence in the 'accessible' region < £ < 1. This correspondence seems to become 
imperfect in the case of non- vanishing subtraction constant As in Ref. [31], we shall argue 
below that this might not be the case, if the £ — > asymptote is given by Regge behavior, as is 
commonly assumed. 

In a real experiment one accesses only a limited phase space region. To pin down CFFs in 
inaccessible regions one can again employ their holomorphic property, but now expressed as a 
family of sum rules. To give an example, we recall that Regge asymptotics of on-shell amplitudes 
yield finite energy sum rules (FESRs), which are widely used in phenomenology. Pragmatically, we 
suppose that off-shell CFFs possess the high-energy asymptotic power-like behavior where 
a(t) is an effective Regge trajectory, which might depend on Q 2 . We introduce now a duality 
parameter v that divides the low- (v < £') and high- (£' < v) energy region in the dispersion 
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integral (fTUl) . Expanding in the vicinity of f = yields FESR0 



S n (v, 



1 

d^- m ~ 2 ^H^---) = ^S m+1 (v,---), with (13) 

— 2 _1 r/ n r i r v ( i i 

For odd (even) m they are valid for signature +1(— 1). We emphasize that the m — — 1 case 
deserves special attention. The resulting FESR is closely related to what is called the inverse 
moment. For real Compton scattering the form of this specific FESR can be found in Ref. [122] . 

The concealed part of the answer to the question 'What GPD information can be revealed 
from experiment?', arises from the analyticity of the amplitude, which ties the real and imaginary 
parts by the single variable dispersion relation ([S]), advocated in Ref. [43]. Plugging the LO result 
([6]) for the imaginary part into the dispersion relation ([8]), it becomes obvious that the real part 
of the amplitude (@J, up to a possible subtraction constant, is also expressed by the outer GPD 
region: 

Ke^(£, d, t, Q 2 ) = PV f dx (— ?— =F -^—] F*(x, V = dx, t, Q 2 ) + C F (d, t, Q 2 ) . (14) 

Jo i + xj 

This 'dispersion relation' tells us that for given \d\ < 1 the variation of the real part with respect 
to £ is governed by the value of the GPD on the trajectory rj = dx. Note that we have introduced 
a new symbol for the subtraction constant fixed to LO: 

C F (d,t,Q 2 ) = C r {d,t,Q 2 ). (15) 

Compared to the LO formula (J4j), in which the whole GPD support contributes to the real 
part, the observation, manifest in Eq. (fl4|) . looks perhaps at the first sight rather surprising. 
It is a consequence of analyticity and the decomposition of the CFFs in terms of GPDs, which 
posses definite properties under Lorentz transformations, clearly visible in the operator product 
expansion approach (OPE) [4"5| 134"] . The consequence is that the GPD in both the central and the 
outer region know about each other and talk to each other via the trajectory t] = dx [83, Hill]- To 
obtain the equality that governs this cross-talk, we combine formula (1141) with the LO amplitude 
(HD and find a family of GPD sum rules (GPDSRs) [431 SS|: 

I dx (— !— =F t^—] [F^(x, n = 0£, t, Q 2 ) - F*(x, V = dx, t, Q 2 )] = C F (d, t, Q 2 ) , (16) 
Jo Vs — x s + x J 



12 Assuming that Q'mJ 7 w ir Pa£,~ a /^{& + 1) is approximately valid for |£| < v one straightforwardly finds 
for a > and n > 0: S n ~ J2 a r(f+i) V a+n ■ ^° obtain the textbook form of S n for on-shcll amplitudes [52], one 
might now rewrite the FESR (JT5]) in terms of the variables v = Q 2 /2M£ and iV = Q 2 /2Mv and take instead of 
the Bjorken limit the on-shell one. To find the more general form, including a < 0, one should follow the lines of 
Ref. |121j . where for n = the a = case has to be treated separately, e.g., as in Ref. [122] . 
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Figure 3: (a) The integration paths in the sum rule (|T6l) along the rj = fix (solid, red) and 
horizontal 77 = (dashed, blue) trajectories are displayed, (b) Visualization of the approximate 
duality relation of a GPD with the s- and t-channel view on the hadronic amplitude. For on-shell 
amplitudes duality between the s- and t-channel view has been conjectured at the end of the 
sixties. Quark-hadron duality was conjectured in particular in DIS, while modelling GPDs in the 
central region by the crossed t-channel contribution was proposed in [SI]. 

parameterized by Here we have dropped the principal value prescription, since the integrand 
can be considered as integrable at x = £. The GPDSR family (Tl6|) is visualized in Fig. [3^. It 
relates the integral along the rj = fix trajectory (solid) in the outer GPD region to a family 
of integrals along the horizontal r\ = i?£ trajectories (dashed), traversing both inner and outer 
regions. The variations of a spectator model GPD on a horizontal (dashed) and the cross-over 
trajectory (solid) are displayed in the left panel of Fig. [TJ 

3.3 Deconvolution and duality 

Let us suppose that the GPD in the outer region and the subtraction constant are known from a 
'measurement'. Finding the GPD in the central region can then be considered as a deconvolution 
problem. Decomposing the integral regions into the central and the outer one, one can view the 
GPDSR family ffl6|) as a Fredholm type integral equation of the first kind, which allows one to 
solve the deconvolution problem. This mapping is an integral transformation that requires the 
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continuation in the region rj > 1. Such an extension is already defined in the DD representatio 

F(x,r ] ,t) = (l-xf [ dy [ * dzS{x - y - r,z)f(y,z,t) . (17) 
Jo J-l+y 

One can verify that the (anti-) symmetrized GPDs within this representation satisfy the GPDSR 
family ( 1161) . The DD representation has been utilized in Ref. [71J to solve the deconvolution 
problem for a simple spectator model, formulated in the light-cone wave function overlap repre- 
sentation, see also Refs. [69], [70] . 

Within the GPDSR family ( |T6l) the deconvolution problem can also be solved by a constructive 
method that explicitly avoids the continuation into the region rj > 1. This method resembles the 
derivation of the standard factorization formula (jlj) from the short distance OPE and the dispersion 
relation ([8]) in the Euclidean region, see Refs. [13, [M] . 

We go now in reverse direction along the lines of Ref. [M]. Expanding the integral kernel in 
the GPDSR family f)16p in powers of l/£ relates integral Mellin moments, 



J dxx^F(x, 7] = i3x, t, jj 2 ) for j > 



of the GPD at the trajectory rj = dx to the usual Mellin moments of the GPD with argument 
r] — 



J dxx j F(x, 7], t, jj 2 ) for j>0. (1£ 



The latter are defined in terms of the expectation values of local operators with spin j + 1. Lorentz 
covariance (or polynomiality condition) enforces 

i+i 

f ] ( V ,t,^) = J2t ) (t,fi 2 )v\ (19) 

n=0 
even 

where ff = for / = {h + e, h, e} (for details see, e.g., O [7]). After reordering, we equate the 
coefficients in front of C,^^ 1 

„1 oo 

/ dxxiF^(x,#x,t,Q 2 ) = V^7j:l(*,Q 2 ), (20) 
Jo 



n=0 
even 



and thus find the coefficients of the polynomial (TT9~j) : 

1 d n 



tf=0 



13 The factor (1 — x) p in our representation (fT?]) is motivated by the form of the positivity constraints [55] and 
appears in the model of Ref. [71]. For the GPD E we have p = 1 and thus a non- vanishing subtraction constant. 
For GPD H + E and the two parity odd GPDs p = and the subtraction constant is absent. 
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where j = 1, 3, 5 • • • and j = 0, 2, 4 • • • for — and + type of GPDs, respectively, and j — n > 
with n even. If a subtraction constant is present in S, the Mellin moments with n = j + 1 appear 
and are evaluated from 

^('■« = b^s^ cw,0>) ~- (22) 



i?=0 



Note that hj +1 \t,Q 2 ) = —ej (t,Q 2 ), and that here, as in the rest of the paper, we follow the 
convention (jSJ), i.e., C = Ce while Cm = —C. 

All possible techniques are now at hand to construct the GPD from the moments (|T9|) . (121]) 
and (I22p : standard inverse Mellin transform, expansion with respect to conformal eigenfunctions 
[59] 160] EH, [42], and introduction of DD as generating function for the Mellin moments [H [72] . 
see, e.g., Ref. |123j . All of them are based or related to the uniqueness of an analytic continuation 
of the even or odd Mellin moments, ensured by the Carlson theorem [124] . Thereby, the Mellin 
moments ej + might either be considered connected with a 'regular part' of the GPD, e.g., as 
in Eq. (fTTj) . or, alternatively, they can be included in the D-term addenda of a GPD [91]. It is 
mostly a matter of taste and convenience (as long as the whole GPD is considered as known), 
which representation one prefers [86]. Only in the case that the spectral function, i.e., SsmJ 7 , 
possesses unexpected analytic properties such as 'invisible' terms, like xS(x), one has to treat the 
Gj terms separately from the DD part. In any case, after convolution of the GPD with the 
hard-scattering part one has the representation for the subtraction constant, cf. Eq. ([22]): 

oo 

C($,t,Q 2 ) = 2^>' +1 e^ +1) (t,Q 2 ) with C($ = 0,t,Q 2 ) = 0, (23) 

3=1 
odd 

where the normalization at d = is the perturbative prediction. This vanishing of the subtraction 
constant can be considered as a consequence of the absence of a local and gauge invariant spin-zero 
operator in the short distance OPE. 

That the GPD in the outer region can be uniquely restored from its knowledge in the central 
one is known and was again clearly spelled out in Ref. [42]. This mapping is (implicitly) done 
in the conformal partial wave expansion approach by taking appropriate boundary conditions for 
the partial waves; see [421 IB!)] for a discussion. Hence, both regions contain the same information 
and so they are dual to each other. 

In anticipation of phenomenological applications we shall now derive GPDSRs that have the 
form of FESRs. Note that GPDs and their sum rules depend on the perturbative approximation 
and scheme conventions that are used in the factorization approach. Analogous to the the deriva- 
tion of FESRs, we introduce the duality parameter v, which divides the integral in the GPDSR 
family ([16]) into the 'low-' and 'high-energy' contribution. By expanding the integral kernel in 
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Eq. ([TBI) around £ = we find the desired GPDSRs: 



dx x 



-m-2 



/ Q \ n j n 

F*(x, Ar, t, Q 2 ) - £ [ ^-^F^(x, V , t, Q 2 ) 



n=0 
even 



5S^ +1 (v,$,t,Q 2 ), (24) 



where m > — 1 is odd (even) for — ( 
obtained from the Taylor coefficients: 



SS: 



d" 



nl 



dp 



dx 



LJQ 



1 



=F 



1 



x 



n\ dif 1 ' ' r]=o 

type GPDs. The 'high-energy' part 5S^(v,$,t,Q 2 ) is 

(25) 



[F T (x, ■ ■ ■ ) - F T (x, tix,---)]- C F (- ■ ■ 

J£=0 

The l.h.s. of the sum rule (1241) relates for x > v the GPD at the trajectory 77 = fix with its 
behavior in the vicinity of 77 = 0, where the strength of the skewness effect is given by its r.h.s., 
i.e., it is embodied in the 'high-energy' constants 8£>*{v ) '&). 

We add that the 'low-energy' and 'high-energy' content might be related by a more elegant, 
however, quite formal way, using analytic (or canonical) regularization [125] . This is useful if one 
deals with GPD models formulated in Mellin space with meromorphic functions. The GPD inherits 
the Regge-like behavior of the CFFs, which we assume, and should not possess any 'invisible' 
terms, e.g., such as x8(x). Under this assumption one can introduce analytic regularization at 
x — 0, indicated by (0) as the lower integration limit, in the sum rule ffTBl) and expand the integral 
kernel in the vicinity of £ = 0. An example is the evaluation of the subtraction constant from the 
'high-energy' limit of Eq. (fT6j) : 



C($,t,Q 2 



lim 



dx 



2x 



[E-(x,#x : t,Q 2 ) - E-(x 7 <&£,t,Q 2 )] 



(26) 



x 2 — £ 2 

It is obvious that interchanging the limit with the integration could render a non-integrable sin- 
gularity at x = 0. Before this interchange, one first has to introduce a regularization, which does 
not spoil the normalization condition C(i9 = 0, t, Q 2 ) = 0. In the absence of 'invisible' terms, such 
as x5(x), we can utilize analytic regularization and write, as in Ref. [34] : 



(27) 



C($,t,Q 2 )= [ dx - \E-(x,$x,t,Q 2 ) -E-(x,0,t,Q' 
J(o) x 



It can easily be verified that this prescription, as well as Eq. (!26j) . is consistent with the GPDSR 
(j24|) for m = — 1. Under the assumption of Regge behavior, an analogous relation should hold for 
the imaginary part of CFFs, too, as argued in Ref. [34] based on knowledge of radiative corrections. 
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4 Phenomenological applications 

4.1 Strategy to access GPDs 

As we saw, the imaginary part of the CFFs in the region x = £ G [0, 1] is a spectral function, 
considered as 'measurable' in our Compton gedanken experiment for < 1. The real part, 
arising from the dispersion relation (jSJ), is then immediately known, too, and it must coincide 
with the 'measured' one in the region x = £ G [0,1]. By means of Eq. (fl2l) . the logic can also be 
reversed: having restricted 'experimental access', i.e., knowing the imaginary and real part of the 
CFFs only in a certain limited interval, one can use the dispersion integral (|S]) (and also Eq. (TTTj) ) 
as a 'local' sum rule, which restricts the shape of the spectral function in the 'non-accessible' 
region. Such a restriction is possible even if one knows the real part only at one given value of 
£. Various forms of sum rules can be derived and we only concentrate on the FESR ones, which 
relate the 'low-energy' and 'high-energy' regions. 'Measuring' the spectral function in dependence 
of for |$| < 1, allows us to solve the deconvolution problem and to reveal the GPD. This can 
be expressed also in terms of GPDSRs in the form of an integral equation ([16]) or in terms of 
FESR-like ones (1241). 

We now apply our results to the DVCS case, in particular, to the dominant CFF 7i: 

0=1, H(U,Q 2 )=n^^ = l,t,Q 2 ), C(£,t,g 2 ) = C(£,0 = l,t,Q 2 ). (28) 

Here we introduce the virtuality of the photon Q 2 « 2Q 2 . Everything that was stated regarding 
the measurement of CFFs remains valid in DVCS kinematics. It is often easier to measure the 
imaginary part (via single spin asymmetries or polarized cross sections) than the real part of the 
amplitude. Hence, one might use the available experimental findings for the real part to reveal 
the imaginary part of the CFF outside of the accessible phase space. 

We introduce two phenomenological parameters 77 and rj to divide the £ interval into the small 
(£ < v)i valence (77 < £ < 77), and large (77 < £) region. We assume that the CFF can be measured 
in the valence region 77 < £ < 77. The dispersion integral (jHJ) might then be used to constrain the 
imaginary part in the large £' region 77 < £' 

^'—^2 [SmH&t, Q 2 ) - S*mH(t,t, Q 2 )] = ttS^(U, Q 2 \%v) + ^ mod (£,t, Q 2 \v) , (29) 
from the experimental knowledge of the region 77 < £' < 77: 

S cxp (£,t, Q%rf) = ReHfot, Q 2 ) ~ \^{t^\ ^fet, Q 2 ) (30) 
+ - f d ^T^-Ti ^m^(r,t,Q 2 )-Sm^(e,t,Q 2 )] , 
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the model estimate for the subtraction constant, and the extrapolation into the small £' region 
i' <ry. 

S mod (£, t, Q\ V ) = C(t, Q 2 ) + - /V-^- [9mW(£', t, Q 2 ) - 3mH(£, t, Q 2 )} . (31) 

^ Jo £ S 

The subtraction constant is known to some extent, see, e.g., Refs. [HI HOOj E3], and one might 
develop some understanding how Regge phenomenology should be employed for off-shell processes. 
In evaluating the amount arising from the very small £ region one can use the DVCS measurements 
of the Hl/ZEUS collaborations [H [T2], [131 HE]- Hence, the model dependent part is therefore 
somewhat under control. The large £ region which one reveals in this way is certainly challenging 
to measure in experiment. As one realizes, even knowledge of the CFF at one value of £, i.e., 
7] = fj = £, leads to a highly non-trivial constraint. This sum rule (|29|) (with £ < fj) can again be 
considered as a family of FESR. To see this, one simply expands the integral kernel with respect 
t° (£/£') 2 • We also note that the real P art of the CFF W can be measured in the crossed channel 
^*^(*) j\rj\r ( or _ > / + /~7) and then obtained by analytic continuation. 

The dispersion integral and crossing symmetry allows for revealing the CFF as a physical 
quantity. The microscopic view on this quantity is given by perturbation theory, providing us, to 
LO accuracy, with a partonic interpretation in terms of the GPD on the cross-over trajectory: 

H~(x, x, t, Q 2 ) = - SmH(C = x, t, Q 2 ) for r]<x<fj. (32) 

7T — 



To that accuracy we write the physical sum rule ([29]) as a GPD one: 
f 1 1r 

J_ dx ^TZ^2 [H~(x,x,t, Q 2 ) - H-{M,t, Q 2 )} ~ C(t, Q 2 ) (33) 
v 2r to 

dx^-^[H-(x,x,t,Q 2 )-H-(C,U,Q 2 )] = S^(U,Q\V,V) for V<i<V- 

It clearly states that in an experiment the GPD can be accessed only on the cross-over trajectory. 

The information which can be revealed away from the cross-over trajectory is governed by the 
GPDSR family fll6p . which we write in the form of an integral equation 

r 1 9r f 1 2r 
PV / dx- -H~(x,rj,t, Q 2 ) = PV / dx- -H~(x, x, t, Q 2 ) - C(t, Q 2 ) . (34) 

Jo T] 2 - X 2 Jo T] 2 - X 2 

The r.h.s. of this integral equation is considered as known and can be equated in an ideal experi- 
ment with the measured real part of the CFF (including the subtraction constant). The l.h.s. is 
nothing more than the real part of the LO result (pjj. In addition, we have the 'boundary condi- 
tion' that H~(x, x, t, Q 2 ) is equal to the imaginary part of the CFF. Knowing the imaginary and 
real part of 7i in an ideal experiment, one would be able to measure the subtraction constant. 
However, the deconvolution problem cannot be solved^ within the family of GPDSRs ( j34t) . 



4 We recall that the scale dependence provides an additional handle on the deconvolution problem. 
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The invisible terms in an ideal DVCS experiment are the zero modes 

2x 



dx- -5H(x,rj,t, Q 2 ) = 0, (35) 

rj 2 — x 2 



in the integral equation ( 1341) with the boundary conditions 

5H(x,rj = x,t, Q 2 ) = 0, lim SH (x, rj, t, Q 2 ) = x^ , lim SH(x, rj, t, Q 2 ) = (1 - xf , (36) 

?J— >0 T)—*l 

where {e, e} > 0. Due to the duality properties, governed by Lorentz covariance, it is clear that 
a priori any given function in the outer or central region which satisfies the boundary conditions 
defines a zero mode. A simple example for fixed Q 2 is 

6H**(x, V ,t) = {r}2 ~ x2) -ff(\x\,t) (37) 
\x\ dx 

with appropriate behavior of f(x,t) at the boundaries x = and x = 1. Hence, the experimental 
data (at least for a fixed scale) can be described by numerous GPD models. In our toy example, 
the zero mode might contribute to the forward case and might be visible for t = in DIS. An 
example of a term visible in DVCS, but invisible in forward kinematics, is given in Ref. [23]. Also, 
for GPD E~ , where phenomenological constraints for the forward limit do not exist, any GPD 
model which fits the data is as good as any other, obtained by a 'gauge' transformation. However, 
the zero mode might contribute to the first moment: 



E~ &t (x, rj, t, Q 2 ) -> E~ Rt (x, 77, t, Q 2 ) + 8E(x, rj, t, Q 2 ) , J dxx8E(x,r],t, Q 2 ) ^ 0. 



Obviously, different GPD models which are able to describe the experimental data can provide 
different values for the desired first moment, which enters the quark angular momentum sum rule 

Although the GPDSR family (|34j) states that the central and the outer GPD region are re- 
lated via the GPD on the cross-over trajectory, some information might be lost in DVCS and/or 
for forward kinematics. For a more geometrical interpretation we refer to the discussion in the 
Introduction and Fig. QXb). Only in the case that for GPDs the holographic principle applies, i.e., 
all information which is contained in the GPD can be obtained from its value on the cross-over 
trajectory, the whole GPD can be revealed from an (ideal) measurement. Nevertheless, the GPD 
sum rules ([3"lj) provide us with valuable information on the 77-dependence of the GPD. 

This can be seen if we decompose the integral into contributions from the central and outer 
regions. Neglecting a possible subtraction constant and assuming that the Regge behavior is 
^-independent, 

P 2x f 1 2x 

/ dx — \H~(x,x,t) - H~(x,r),t)] = / dx — [H~(x, x, t) - H~(x, 77, t)\ , (39) 

J rf - x 2 L J ^ x Tj J 

2(3 



we find that in the limit 77 — > the l.h.s. vanishes and the r.h.s. must therefore satisfy an inverse 
moment sum rule. Taking now the limit 77 — ^ 1 and using the fact that the GPD is vanishing at 
the end point x — 1, we find another sum rule. These sum rules read 

[ dx - [H-(x,x,t) -H~{x,0,t)] = 0, [ dx X 2 [H~{x, x, t) - H~{x, 1, t)] = . (40) 
Jo x Jo 1 — x 

In the left sum rule (i.e., 77 — 0) the GPD on the cross-over trajectory is dictated by the physical 
content in the outer region, while in the right sum rule (i.e., 77 — > 1) it is governed by the one in 
the central region. This suggests that for small 77 we shall naturally have the partonic s-channel 
point of view, while for large 77 the t-channel point of view in terms of meson-like exchanges is 
more appropriate [S3]; see left and right sketch in the lower part of the Fig. [3]^b), respectively. To 
match both regions one might sum up the t-channel exchanges, i.e., reggeize them. If one would 
like to incorporate phenomenological knowledge of hadronic physics in GPD modelling, one might 
suppose that for small x and 77 one understands better the hadronic t-channel view in terms of 
mesonic Regge exchanges [33], while for large 77 and x the s-channel partial wave expansion in 
terms of few baryonic resonances is more useful [127] ; see right and left sketch in the upper part of 
the Fig. E^b), respectively. Hence, GPD modelling from phenomenological input should be related 
as indicated across the four sketches in the Fig. [3](b) rather then vertically. Finally, we note that 
the left equation in Eq. (|4*0|) is given in its more general form in Eq. (149|) . 

Let us also discuss the problem of relating the GPD at the cross-over trajectory with the zero 
skewness GPD from the point of view of Mellin moments. The subtraction constant that enters 
the DVCS process is, according to Eq. (j23l) . given by 



C(t,Q 2 ) = -2^^ +1) (t,Q 2 ). (41) 
j'=i 

odd 

A model dependent estimate of the constant for t — was provided in xQSM [73]. The xQSM 
model evaluation of the first term in the series (|4ip . given in Ref. [100] . is compatible with lattice 
measurements [83]. This allows one to estimate the size of the constant C(t, Q 2 ), to judge on the 
convergence of the series (l41j) . and to estimate its t-dependence. According to Eq. (|20|) . the Mellin 
moments of the GPD at the cross-over trajectory are given by the series 

„! 00 

/ dx x?H-{x, x, t, Q 2 ) = hfl(t, Q 2 ) . (42) 

^0 n 



n=0 
even 



The first term with n = in this series is the j th Mellin moment of the GPD at rj = 0. 

Instead of using the number of total derivatives n as a label for the degrees of freedom, one 
might prefer the angular momentum J of the t-channel SO (3) partial wave expansion, which is 
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a physical degree of freedom. If we apply crossing symmetry to the GPDSR family ([Mil , for 
details see, e.g., Ref. [12], and set r] 2 1/ cos 2 0, where # is the scattering angle in the t-channel 
center-of-mass frame, we obtain the LO approximation of the corresponding 77* — > pp amplitude 
in terms of the generalized distribution amplitudes H^~(z, cos 9, s, Q 2 ) pi [85]: 



fdz -^- 2 H^- (z, cos 9, s, Q 2 ) = fdx 2xC f 6 2( . H-{x, x, t, Q 2 ) - C(t, Q 2 ] 
Jq 1 — z J 1 — x cos o 



t-*a 



(43) 



In the CFF combination Ji + t£/4M^ Legendre polynomials Pj (cos 9) with angular momentum 
J = 0,2,4, ••• are the proper partial waves [6]. The corresponding partial wave amplitudes 
dj(s, Q 2 ) can be obtained from the orthogonality relation for Legendre polynomials. In other 
words, one multiplies the l.h.s. or r.h.s. of, e.g., Eq. ([43]) with the partial waves and integrates over 
Z = cos 9. It is then clear from Eq. (1430 that the subtraction constant C(s, Q 2 ) contributes only 
to the J = partial wave amplitude. For J > 2 the r.h.s. of (143j) tells us that the partial wave 
amplitudes can be obtained from the Froissart-Gribov projection by analytic continuation in t: 



a As, Q ) = 2 / dx 



x 2 



t 



H~{x, x, t, Q 2 ) + Tjjz-E-ix, x, t, Q 



±M 2 N 



(44) 



where Qj{z) are the Legendre functions of the second kind. The representation ( l44l) is suited for 
the analytic continuation in J within a(t) < KeJ, where a(t) is the leading Regge trajectory. The 
partial wave amplitudes are now, via Eq. ([44"]) . represented in terms of partonic quantities, namely, 
by the GPD on its cross-over trajectory. We might also express them as a series of common Mellin 
moments, labelled by the spin and the number of total derivatives, or by the moments with definite 
conformal spin [34] , 

Obviously, the perturbative approach 'predicts' the partial wave amplitudes in terms of GPDs. 
Since the GPDs are unknown, it is more appropriate to employ the reverse logic. As far as the 
partial wave amplitudes are known they provide a strong constraint on the GPDs. A power-like 
parameterization of the GPDs at its cross-over trajectory at small x, i.e., x~ a ^\ corresponds 
to a pole at J = a(t) in the SO(3) partial wave amplitudes and makes direct contact to Regge 
phenomenology, extended to off-shell kinematics. For integral values of J > 2 we can read Eq. ([4*4"]) 
as a sum rule. The case J = deserves special attention. It can be unproblematically evaluated 
from Eq. (1431) . Naively setting J = in Eq. (j44p provides a divergent integral as the natural 
outcome of our standard derivation. We should require that both J = results are equal, which 
finally would provide a definition of the inverse moment for the GPD at its cross-over trajectory. 
Such an inverse moment is then simply given by the subtraction constant C(t, Q 2 ). 

The GPDSRs, such as in Eq. ( 134]) . might be applied in various manners and we have presented 
here only a few illustrative examples. Finally, we would like to point out that they also make 
contact to the work of Refs. [H HQ] ST]. Reducing the GPDSRs dSJ to the DVCS case, i.e., 
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$ = 1, we find in the limit v — > 0: 



1 



^TlH - 1 yi YFL 2 ,-772 

2 i^(*, x, t, g 2 ) - £ F^, ^ t, g 2 

n=0 ' ' 

even 



rj=0 



\im6S^ +1 (v,t,Q 2 ), (45) 



u->0 



where S^ +1 (v,t, Q 2 ) = S^ +1 (v,$ = l,t,Q 2 ) can be read off from Eq. (125]) . If we now view all 
terms as analytic functions of m, we can invert this specific GPDSR (145 p . by an inverse Mellin 
transform to express the GPD at its cross-over trajectory in terms of 'low-', i.e., second term in 
the square brackets of Eq. (1451) . and 'high-energy' constants, i.e., the terms of This can 

then also be converted in a series of forward-like parton distributions q u (x,t, Q 2 ) weighted with 
x u . The expansion of the GPD on its cross-over trajectory in this forward-like parton distributions 
can be written, e.g., for H~, aJ^I 

oo 

H~(x, x, t, Q 2 ) = ^Tx v q„(x, t, Q 2 ) with q (x,t, Q 2 ) = H~(x,r] = 0,t, Q 2 ) . (46) 

even 

Here the forward-like parton functions know about both the 'low-' and 'high-energy' part and 
should be viewed as generalized functions in the mathematical sense. The polynomiality condition 
translates into 

/ dxx™' 1 q n (x,t, Q 2 ) = for 2 < m< n , (47) 
Jo 

where m and n are even. The decomposition of these forward-like parton functions in forward-like 
conformal or 'dual' parton functions has been discussed in the small rj expansion in Ref. [126] . 
Generally, it should be given by an integral transformation, cf. Refs. [601 EH EH SOI SI]- We 
emphasize that these forward-like functions, e.g., given by explicit expressions in Ref. [41], can be 
interpreted in terms of GPDSRs (1451) . The DVCS amplitude to LO in the 'dual' parameterization 
is then finally represented as a sum of functions that are associated with the 'low-' and 'high- 
energy' content, where the duality scale v itself is taken to be zero; see Ref. [41]. This step 
generates the 5-functions (including derivatives) which appear in the 'dual' parameterization. We 
recall that in the series (]46p all forward-like parton distributions contribute to the small x region 
and that at least the first two terms are needed to describe DVCS to LO accuracy in the small 
XBj region. 



15 One might find this result directly from the 'Taylor expansion' of H~(x,r],t, Q 2 ) in the vicinity of r/ = with 
x = 7/. However, the expansion coefficients are generalized functions, whose property (47) ensures polynomiality. 
These forward-like parton distributions can be straightforwardly evaluated from DDs and the expansion (TIT)]) can 
be studied, e.g., within a RDDA. 
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4.2 Revealing GPDs from experimental measurements 

We address now the question: "What do we learn from a DVCS measurement?". To have a 
partonic interpretation one is mostly interested in the GPD at t] = 0. The extrapolation of 
F T (x,x,t, Q 2 ) to F T (x, 0,t, <2 2 ) is a model dependent procedure, which can be parameterized by 
the skewness function S(x,t, Q 2 \F :f ) defined by 

F T (x, x, t, Q 2 ) = [l + S(x, t, Q 2 |F T )j F T (x, rj = 0, t, Q 2 ) . (48) 

Note that we included a t-dependence in the skewness function, since it is expected that this 
dependence at large x is rather different from that of the GPD at 77 = 0. The GPDSR family (fT6|) 
provides us in the limit £ — > with a link between the measured GPD at 77 = x and 77 = and 
therefore a constraint for the skewness function. Since the limit might be intricate, we first have 
to regularize: 



/ dx -S{x,t, Q 2 \F-)F-(x, V = 0,t, Q 2 ) = \c F {t, Q 2 " 
J(o) x 2 



(49) 



One might consider the F^(x,x,t, Q 2 ) as new functions, which are in an unknown manner 
connected to the GPD at 77 = 0. As emphasized above, the analytic regularization procedure 
might be reinterpreted in terms of partonic FESR. The constraint ()49p on the skewness effect is 
rather weak; only the inverse moment is constrained, and so there are infinite degrees of freedom 
left, e.g., in the series (|4ip . We might parameterize the skewness effect by relative deviation 
factors: 

s(t W _ y fjtM _ tidxxiS(x,t,f, 2 \FT) F T(x, v = 0,t,/i 2 ) 

hi /i(*^ 2 ) J 1 dxx^(x,77 = 0,t,/i 2 ) ' 1 ' 

even 

with, as before, j being odd (even) for the — (+) superscripts. The Mellin moments of CFFs are 
then given as a product, 

[ dt¥%mF(U,Q 2 )=K /Wei 7lT (^ = e,e,t,Q 2 )=vr/,(t,Q 2 )[l + ^(t,Q 2 |F^)] , (51) 
Jo Jo 

of the desired Mellin moment fj(t, Q 2 ), surviving the forward limit, times the skewness effect. 

For t = the Mellin moments hj(t = 0, Q 2 ) [hj(t = 0, Q 2 )} are those of the PDFs and are 
constrained by phenomenological input. The lowest Mellin moments fj=o{t, Q 2 ) are partonic form 
factors, related to observables. The deviation factors are built by a series (1501 and the present 

state of the art in lattice measurements, to the best of our knowledge, allows for the evaluation of 

(2) 

the first term in this series for j = 0, i.e., the quantity h 2 . As long as the forward limit is known 
(to some extent) a DVCS measurement of, e.g., TL allows us to pin down the skewness effect for 
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this quantity. In the case one wants to extract, for instance, the quark angular momentum in an 
ideal DVCS experiment via the anomalous gravitomagnetic moment [5], 

ti\M ) - ei(t — u, y j — ^ 1 + ^(^ = 0,22^-) > ^ 

one certainly needs to know the skewness effect (at t = 0), which might be different from that 
appearing in H. 

The conclusion is obvious. One must use additional information to reveal the quantities one 
likes to pin down. Lattice measurements can give only a hint on the size of the deviation factors, 
and this we consider very valuable. Model evaluations, e.g., within xQSM, can at least give 
estimates for the size of deviation factors. As in any of such model calculation, the matching of 
collective with current quark degrees of freedom induces uncertainties. Nevertheless, it is very 
important to quantify these deviation factors on a generic level. Further experimental constraints 
on them, available in the small xbj region, arise from the logarithmical scaling violation, governed 
by evolution. A discussion of revealing the quark angular momentum from the present DVCS data 
should be given elsewhere. We would only like to note here that GPD model dependent constraints 
on the quark angular momentum require understanding of the skewness effect. A minimal basic 
requirement for such a constraint is that the utilized GPD model describes all available DVCS 
data, including the presently challenging data on unpolarized cross section. 

According to the above discussion, one strategy for extracting GPD information to LO accuracy 
might consist of the following steps: (suppose that we consider H or H) 

• model the t-dependency of the GPD at 77 = 0, 

• parameterize the skewness function, defined in Eq. ([48]) . and 

• fit the parameters to the measured observables (CFFs) by means of Eqs. (1321) and (1331) . 

Alternatively, in the third step one might use the sum rule (1291) to first pin down the imaginary 
part in the whole region from knowledge of the real part and then reveal from Eq. fl32l) the GPD 
at the cross-over trajectory. 

Certainly, in the first step it is crucial to have a realistic model for the interplay of the x- and 
t-dependence of the rj = GPD. Proposals for the functional form of the GPD at 77 = are given 
in Ref. [1281 11291 |9~T] . We would like to suggest another empirical recipe that allows one to take 
advantage of theoretical/phenomenological knowledge on the generic level. Moreover, it has a 
power-like behavior at large —t and can in principle also be used for positive values of t. Inspired 
by the simplicity of the Veneziano model [130] for on-shell amplitudes, we found it simple to set 
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Figure 4: Panel (a) shows various models for the dominant u-quark GPD H at rj — x versus the 
Bjorken-like variable X = 2x/(l + x) for fixed t = —0.25 GeV 2 . In panel (b) we confront these 
models with JLab/Hall A measurements of the observable QmC^,, defined in Eq. (1541) . for DVCS 
off unpolarized proton [TIj for a^Bj = 0.36 and Q 2 = 1.9 GeV 2 . 



up GPD models at r/ = in Mellin space by taking ratios of T functions, e.g., for the valence 
c?-quark contribution: 

. ^ = r(i + j - q(t))r(i +j + p- a)r(2 -a + p) 

j T(l-a)T(l+] + P-a(t))T(2 + 3-a + f3)' { ' 

The generic values of the parameters arise from counting rules (j3 = 3, P = 2), Regge phenomenol- 
ogy (a(t) = 1/2 + t/GeV 2 ), and for u-quarks we employ in addition SU(6) symmetry arguments, 
i.e., h J val / h^ val = 5 for j — > oo. These parameters are then understood as effective ones and 
are slightly tuned to phenomenological PDF parameterizations and form factor measurements. 
Lattice results [83] can be quantitatively reproduced within our ansatz if we adopt the found 
parameterization of the Regge trajectory a(t) to the large meson mass lattice worldo The guid- 
ance for the parameterization of the skewness function S(x,t, Q 2 = 2GeV 2 |if~) comes from the 
functional form of the solution for the positivity constraints [50], the expected large x behavior 
for GPDs at rj = x [131] . and model results [71] . The skewness function satisfies the sum rule 
fj4^|) . where the t-dependence of the subtraction constant has been taken from the first coefficient 
of the D-term from the xQSM result [100J. The small x behavior of the GPD H is taken from 
our LO fits to the Hl/ZEUS data [88] . The modelling will be described in detail somewhere else. 

In Fig. H^a) we show various versions of the model for the dominant w-quark GPD H . Only 
the solid curve is rather similar to the model versions of Ref . [27] ; the two other ones belong to a 



16 We do not share the wide-spread belief, e.g., in Ref. [55], that present lattice data should be used as quantitative 
GPD constraints. 
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new class of GPD models. As argued in Sect. 12.2.31 they are enhanced in the large X = 2x/ (1 + x) 
region. Note that such a functional form resembles the shape of the energy spectrum function 
of an exclusive process, if 1/X is viewed as a rescaled energy (perhaps one would even expect a 
resonance structure before the 'continuum' starts). The w-quark GPD H is the dominant one in 
DVCS off an unpolarized proton target and it might be revealed from the azimuthal dependence 
of the interference term. The quantity which determines the magnitude of the dominant first 
harmonic in the interference term [27] , 

cl p = F^H&t, Q 2 ) + £ [*i(f) + F 2 (t)]H(U, Q 2 ) - j^pMt)£(U, Q 2 K , (54) 

is shown in Fig. 11(b) for JLab/Hall A kinematics: 

x Bj = 0.36 , Q 2 = 1.9 GeV 2 , and 0.15 GeV 2 < -t < 0.35 GeV 2 . 

The relative contribution of the terms involving CFFs 7i and £ we estimate to be on the level 
of 10%-40% and up to ±10%, respectively, compared to the dominant H, term. Instead of ex- 
ploiting the sum rule family fl29l) . we confronted our GPD model with experimental data in the 
old fashioned way. We varied the unknown contributions and compensated for this by tuning the 
skewness function for H. We emphasize that any curve in Fig. IH(a) that crosses one of the shown 
ones at X = 0.36 defines a priori a GPD model, which is not ruled out by the JLab/Hall A data 
in panel (b). The GPD models are then selected by confronting them with the measurement for 
the real part for the observable (1541) . We found that the extracted values for this quantity are 
contaminated by DVCS squared terms and so there is a further contribution from the GPD E. 
Taking all this into account we had no problems to tune the model parameters to the measured 
values for t = —0.25 GeV 2 . However, we were only able to get a qualitative agreement with the 
measured t-dependence of this quantity. Whether this is an artifact of our own GPD model re- 
strictions or is due to the power suppressed terms in the twist-three formalism of Ref. [27] remains 
for the moment an open question. This is in our opinion more a technicality than an issue, which, 
of course, must be investigated. 

5 Outlook 

We have considered CFFs as holomorphic functions and discussed the phenomenological appli- 
cation of this property. To LO accuracy and for fixed photon virtuality one can in a DVCS 
measurement access the GPD only on its cross-over trajectory t] = x. Combining dispersion re- 
lation and perturbation theory allows for a derivation of the partonic analogue of FESR, which 
relate the small and large x behavior of the GPD. Altogether, this sets up a new strategy for re- 
vealing GPDs from experimental DVCS measurements, utilizing well-known techniques from the 
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phenomenology of on-shell hadronic processes. The proposed alternative LO GPD approach has 
numerous advantages with respect to the present ad hoc GPD model approach. The LO GPDs 
can be considered as 'super' GPDs that contain all perturbative and non-perturbative corrections. 
If the assumed holomorphic properties hold true, the approach by definition cannot fail to describe 
DVCS data for fixed photon virtuality. All discussions about GPD models/ansatze and erroneous 
conclusions coming from their success or failure are avoided from the very beginning. The goal of 
this concept is to reveal GPDs at their cross-over trajectory rj — x only from DVCS data and to 
obtain a generic understanding of the skewness effect. 

The partonic interpretation of the GPD to LO accuracy has been given here only from the prin- 
cipal point of view, where the appropriate choice of scaling variables was not discussed. This rather 
important problem is related to the minimization of not well-known higher twist contributions, by 
appropriate variable definitions. We recall that a partonic interpretation of the 'measured' GPD 
at the cross-over trajectory has already been pointed out in Ref. |132j . We would like to specialize 
this interpretation. In the infinite momentum frame one might label the proton state by the center 
of plus-momentum, given as the expectation value of the two generators for transversal Galilean 
boost [133j . The emission of a parton with momentum fraction X = 2x/ (1 + x) leads to a change 
of the center of plus-momentum 6™ of the initial nucleon. Since for 77 = x the outgoing spectator 
system absorbs a quark with zero longitudinal momentum fraction, its center of plus-momentum 
will not change and it coincides with the one of the final state nucleon b s f ec = b^_ ut . The difference 
Sb± = b" 1 — b s f cc can now be viewed as the transversal distance of the struck parton from the 
collective spectator system. One might now adopt the GPD wave function overlap interpretation 
in the impact parameter space [134] . which tells that the Fourier conjugate variable of Sb^ is 
simply the transverse momentum transfer Aj_, directly observed in experiment. In other words, 
following the procedure of Ref. [132j gives us access to the transversal separation distance of the 
struck quark and collective spectator system^ A probabilistic GPD interpretation as for the zero 
skewness case, however, does not hold here anymore. 

Also, the skewness effect can be addressed from the experimental side. In DVCS observables 
for an unpolarized nucleon target one might first neglect the non-dominant CFFs H., S and S 
as well as twist-three contributions. This seems to be justified to a certain extent, at least for 
HERMES kinematics. For the dominant quantity H the corresponding GPDs in the zero skewness 
case are well constrained by PDFs and form factors. Guidance from model calculations and lattice 
data might be employed with care, too, to set up field theory inspired GPD models for the zero 
skewness case. The skewness function, which relates the GPD for zero skewness to the one on the 

17 This interpretation has been also presented in Ref. [135], including considerations on the i-behavior of GPDs 
at large -q = x as it arises also from powerdike wave function models [Ml EH EI] ■ 
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cross-over trajectory is constrained by GPD sum rules (GPDSR). It can then directly be fitted 
to experimental data. The limitation of this pragmatical approach for the collider experiments 
HI and ZEUS is also obvious: the perturbative prediction of the scale change, which gives us an 
additional handle on the skewness effect, cannot be implemented. Since the functional form of 
CFFs and GPDs drastically simplify in this kinematics, experimental data can be analyzed in the 
traditional manner within a flexible GPD parameterization in the small x region. This can be 
done in any GPD representation; however, it has up to now only been applied in the conformal 
approach. 

We expect that a global LO analysis of present DVCS data can already provide a generic un- 
derstanding of the GPD at the cross-over trajectory, in particular for the dominant GPD H. Once 
this understanding is reached one can address the deconvolution problem, i.e., to find a realistic 
parameterization of GPD models. This is particularly important if in dedicated experiments the 
CFFs £ and £ are measured and one wants to access the GPDs E and E, where the former one 
gives access to the quark angular momentum. Realistic and flexibly parameterized GPD models 
should then open the road to the inclusion of radiative corrections and the Wandzura-Wilczek 
approximation of twist-three GPDs, see, e.g., Ref. |136[ 1137] . or to address quark-gluon-quark 
GPDs |138j . If experimental data for virtual electroproduction of mesons are analyzed within the 
GPD formalism, it might be inevitable to go to next-to-leading order accuracy. The formalism 
has been worked out for the twist-two sector (TDTJ, [TDSJ EU 1135111301 ■ 

In this paper it has clearly been spelled out that GPDs have dual properties. If CFFs possess 
Regge behavior, this duality relation is perfect and the GPDs contain the same information in 
the central and outer region. This duality property is implemented in a GPDSR family and can 
be expressed in terms of GPD moments, too. Unfortunately, these GPDSRs are not sufficient to 
find a unique GPD model from (ideal) experimental DVCS measurements. There are two possible 
ways to get an handle on the remaining degrees of freedoms. 

The first way relies on the common QCD tools: phenomenological knowledge, counting rules, 
model calculations, and lattice simulations. Certainly, lattice simulations can only provide a flash 
on the skewness effect, while model calculations suffer from the matching problem of partonic and 
collective degrees of freedom. A central theoretical task here is to understand the small x behavior 
of GPDs. It seems to be rather promising to study the t-channel ladder, e.g., to proceed along the 
lines of Refs. jS] or |117] . We emphasize that phenomenological knowledge can be implemented 
from the i-channel point of view for small x and also from the s-channel point of view for large 
x, where both regions are tied together by GPDSR. This suggests to look for dual GPD models 
in which the meson-baryon duality conjecture, employed in the phenomenology of hadronic on- 
shell processes, is implemented. Crossing relations via the Froissart-Gribov projection provide an 
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additional phenomenological handle. 

The second way to address the problem is to develop a theoretical understanding of whether 
there is some holographic principle that tells us how the information about the GPD on the 
cross-over trajectory is related to its values in the central (or outer) region, which includes the 
rj = case. This question arises naturally if one relates the partonic GPD interpretation to the 
t- and s-channel hadronic point of view. We recall that the conjecture of duality originated in 
experimental findings about hadronic on-shell processes and was extensively studied within FESR; 
see, e.g., Ref. [52]. Thereby, the SO (3) partial wave expansions in both the t- and s-channel were 
utilized in the context of Regge phenomenology, which provides an empirical description of the 
strong interaction phenomena. It was also realized that duality phenomena could be partially 
explained within the quark model. In dual hadronic resonance models the underlying symmetry 
of duality is the SO(2,l) group. Finally, the quark model was gauged providing us QCD, and string 
theory was born out of dual hadronic resonance models, i.e., out of the group SO(2,l). Certainly, 
in QCD phenomenology quark-hadron duality, in particular duality between the deep-inelastic 
and the resonance kinematical regime (Bloom-Gilman duality) is an important phenomenological 
concept; for a review see Ref. [141j . Also conformal symmetry shows up in various QCD corners 
|142[ I143[ I144[ 11451 11461 [58] and integrability phenomena in the multi-color limit were discovered in 
QCD [HH [HBl [HH1 [inni HSU and in supersymmetric Yang-Mill theories; see, e.g., Refs. [T52l[T53l 
and references therein. These integrability phenomena play a key role in getting a handle on the 
AdS/CFT conjecture [531 EH ES]. Initiated by this conjecture, model dependent approaches to 
a supposed AdS/QCD duality were developed [1541 1155[ 1156] . see also Refs. [1571 1158] . and its 
application to high-energy physics is discussed, e.g., in Ref. [159j . From this more historical line, 
or, better to say, loop, we believe that it is rather worth it to study GPDs from the AdS/QCD 
point of view. Certainly, it is not expected that there is an exact answer for (scheme dependent) 
partonic quantities; however, one might perhaps find guidelines for the parameterization of GPDs 
with respect to their skewness dependence. 
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